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Abstract. A recent unambiguous detection of X-ray rotational modulation of the supersaturated star VXR45 (P 
— 0.223 days) has shown that its corona has discrete dark and bright X-ray regions. We suggest that due to the 
rapid rotation, the X-ray emitting corona has been centrifugally stripped away, creating open field regions that 
are dark in X-rays. This leads naturally both to a significant rotational modulation in X-rays but also to the lower 
X-ray luminosity of supersaturated stars compared to those rotating more slowly. To demonstrate the effect, we 
take as an example a more slowly rotating star for which surface magnetograms are available. We extrapolate 
the potential coronal magnetic field based on these magnetograms and determine for a hydrostatic, isothermal 
atmosphere the structure of the density and of the optically-thin X-ray emission. We show that if the rotation rate 
of this star were increased, the magnitude of the X-ray luminosity would decrease while its rotational modulation 
would increase in a way that is consistent with the recent observations of VXR45. 
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While the nature of the X-ray emission from the Sun is 
relatively well understood, the way in which that emis- 
sion changes with increasing rotation rate is less clear. 
Observations suggest that, compared to the Sun, stars 
of increasing rotation rate show a rise in their X-ray 
emission that reaches a maximum of about Lx/Lboi = 
10~^ at rotation rates of about vsini = 15 — 20kms~"'^ 
HVilhu 1984|l . Beyond this rotation rate is the "saturated" 
regime where the X-ray luminosity is independent of ro- 
tation rate. This behaviour persists until rotation rates 
of about wsinz > lOOkms"^, where the X-ray luminos- 
ity begins to decrease again. This regime is referred to as 
"supersaturated" (IProsser et al. 19 96 Randich 1998|. 

There are several possible explanations for saturation 
of the X-ray emission. It may be that with increasing ro- 
tation rate, the dynamo process itself saturates due to the 
back-reaction of the field on the plasma. Alternatively the 
increasing dynamo activity may lead to a complete cov- 
erage of the stellar surface in active regions such that no 
further increase in X-ray emission is possible IjVilhu 1984|l . 
A rise in coronal temperatures above the detection limit 
might also produce an apparent saturation. 

A detection of rotational modulation of X-ray emis- 
sion in a saturated star could potentially address some of 
these questions. A lack of any modulation could be con- 
sistent with a corona that is densely packed with X-ray 
emitting loops, while the presence of modulation might 
provide some clues as to which magnetic structures are 



dominating the emission. Rotation modulation has, how- 
ever, proved very difficult to detect, given the high levels 
of X-ray variability in active stars. 

There have been some indications from observations 
of both binary and single stars that rotational mod- 
ulation may be greater for lower-temperature plasma 
which has a smaller pressure scale height and so may be 
confined closer to the stellar surface IjWhite et al. 1990l 
IGiidel et al. 1995|l . Lack of rotational modulation may of 
course be due to an extended, densely-packed corona, but 
it can also be due to most of the emission coming from high 
latitudes where it is not eclipsed rSiarkowski et al. 19961 
ISolanki et al. 19971 IJeffries 1998. .St§pieh et al. 2001| ). A 
very clear example of this is the case of AB Dor (P=0.514 
days). A long term ROSAT study HKiirster et al. 1997|) 
showed a very low rotational modulation of only 5 — 
13%, while observations of two flares with BeppoSAX 
fMaggio et al. 20001 ) showed no rotational modulation of 
the X-ray emission over the decay phase of the flares, al- 
though they lasted for more than one rotation period. 
Modelling of the flare decay indicated that the flaring 
loops were small, with a maximum height of only O.SR^,. 
This implies that the flaring regions must have been lo- 
cated at latitudes above 60° (the stellar inclination) where 
they were never eclipsed. 

Zeeman-Doppler images of AB Dor 
(IDonati fc Colher Cameron 19971 IDonati et al. 1999|l 

certainly show flux at the kilogauss level at these lati- 
tudes. By extrapolating the coronal magnetic field of these 
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magnetograms and determining the corresponding X-ray 
emission from an isothermal, hydrostatic atmosphere, 
Jardine et al. (2002) showed that the surface magnetic 



maps would lead naturally to most of the emission coming 
from the high latitude regions. This gave a low rotational 
modulation coupled with the observed high density and 
emission measure. 

In the case of the supersaturated star VXR45, how- 
ever, the X-ray emission is very different. This star 
is very rapidly rotating member of IC 2391 with a 
spectral type of dG9 and photometric period of 0.223 
days l|Patten fc Simon 1996|l and a vsini > 200kms~^ 
IjStauffer et al. 1997|l . Its value of log(Lx/iboi) is between 
—3.60 and —3.62 and so it is clearly in the supersatu- 
rated regime. The XMM-Newton observations reported 
by Marino et al. (2003) show an unambiguous rotational 
modulation of about 30%. Given the broad energy band- 
pass of XMM-Newton, [Marino et al. (2003)| rule out the 
possibility that the temperature range of the emission is 
responsible for the supersaturation and claim that the re- 
sult shows clearly that the X-ray emission is coming from 
discrete structures that are rotationally self-eclipsed. 

It is possible that for these very rapid rotators the 
dynamo process has changed its nature from a satu- 
rated state to one where it produces a more patchy 
coverage of flux at the stellar surface. Doppler images 
of VXR45 do not appear to support this suggestion 
IjMarsden et al. 2003|l . but even if this were the case, the 
effect of coronal stripping HJardine fc Unruh 1999|l could 
mask the onset of such a change, or disguise it completely. 
At high rotation rates the co-rotation radius (where cen- 
trifugal forces balance gravity) moves inside the X-ray 
emitting corona. The rise of gas pressure in the summits of 
magnetic loops then breaks open these loops to form open 
field regions that are dark in X-rays. This reduction in the 
emitting volume initially balances the rise in the density of 
the corona to give a saturation of the X-ray emission, but 
eventually enough of the coronal volume has been forced 
open that the X-ray emission falls with rotation rate. As a 
result, at the highest rotation rates much of the corona is 
filled with open field and so there should be a significant 
rotational modulation in X-rays. 

The purpose of this paper is to demonstrate that a 
simple increase in the stellar rotation rate is sufficient to 
produce both a rise in rotational modulation and a drop 
in the emission measure consistent with that seen in both 
AB Dor and VXR45. 

2. Modelling the coronal emission 

We do not have a map of the surface magnetic field of 
VXR45, although Doppler images ( Marsden et al. 2003|l 
show a distribution of dark spots essentially similar to 
that found on the more slowly rotating star AB Dor, 
with much of the surface, especially the high latitude re- 
gions, showing large spots. Given that we do not know 
in detail how the surface flux distribution should vary 
with rotation rate for these saturated and supersatu- 




LONilTUDE IN DECREES 

Fig. 1. A map of the surface radial magnetic field of AB 
Dor. White represents — 800G and black represents 800G. 
Since AB Dor is inclined at 60° to the observer, Zeeman- 
Doppler images provide only limited information in the 
lower hemisphere. In order to compensate for this, we have 
generated this combined surface map, with the 1995 map 
in the upper hemisphere and the 1996 map in the lower 
hemisphere. 

rated stars, we choose to take one sample surface flux 
distribution and keep that independent of rotation rate. 
To this end we have used a magnetogram for AB Dor 
to provide a coronal field structure with the degree 
of complexity that is implied by Zeeman-Doppler maps 
IjDonati fc Colfier Cameron 1997IIDonati et al. 1999|l . We 
note that due to the stellar inclination of AB Dor, only one 
hemisphere can be observed and so to allow for the fiux 
from the hidden hemisphere we artificially add in a surface 
map from another year (see Fig^. Using this composite 
surface map as a boundary condition for the magnetic 
field, we can extrapolate the coronal field by assuming it 
to be potential. 

The method of extrapolating the coronal field has been 
described 
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not be repeated in detail here. We 
source surface method pioneered by 



Altschuler & Newkirk, Jr. (1969) and a code originally de 



veloped by van Ballegooijen, Cartledge & Priest (1998) 
Briefly, we write the magnetic field B in terms of a fiux 
function "if such that B = — and the condition that 
the field is potential (V x i? = 0) is satisfied automati- 
cally. The condition that the field is divergence-free then 
reduces to Laplace's equation V^^* = with solution in 
spherical co-ordinates (r, 9, (f) 
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(1) 



where the associated Legendre functions are denoted by 
Pim- The coefficients aim and bim are determined by im- 
posing the radial field at the surface from the Zeeman- 
Doppler maps and by assuming that at some height Rg 
above the surface the field becomes radial and hence 
Be{Rs) — 0. This second condition models the effect of 
the plasma pressure in the corona pulling open field lines 



M. Jardine: Coronal stripping in supersaturated stars 



3 




Fig. 3. Shown is the emission measure based on the mag- 
netogram shown in FigQ] Results are shown for two as- 
sumed stellar inclinations, 90° (solid) and 60° (dashed). 

to form a stellar wind. We set the source surface to be 
at the Keplerian co-rotation radius such that for a stellar 
rotation rate uj, = (GM/w2)i/3. Thus, = 2.71i?^ 
for AB Dor and 1.55i?^ for VXR45. 

In order to calculate the X-ray emission, we must de- 
termine the coronal density structure. As a first step, we 
calculate the pressure structure of the corona assuming it 
to be isothermal and in hydrostatic equilibrium. Hence 

the pressure at any point is p = pgc™^'^"^ ""'^^ where 
gs = (g.B)/|B| is the component of gravity (allowing for 
rotation) along the field and 



g{r, 0) = {-GM^,/r^ + oj^r sin^ 9, u^r sin 9 cos ( 



(2) 



At the loop footpoints we scale the plasma pressure po to 
the magnetic pressure such that Po{9, (p) — KBq{9, <p). We 
set K = 10~^ for all models in this paper, giving a pres- 
sure of 120 dyne cm^"^ for AB Dor at lO^K based on an 
emission-measure weighted density of 4.4 x lO^'^cm^'^. As 



described in Jardine et al. (2002) scaling K simply scales 
the mean density and the emission measure. The plasma 



Fig. 4. Shown is the rotational modulation of the X-ray 
emission measure based on the magnetogram shown in Fig 
n Results are shown for two assumed stellar inclinations, 
90° (solid) and 60° (dashed). 



pressure within any volume element of the corona is set 
to zero if the field line through that volume element is 
open, or if the gas pressure exceeds the magnetic pres- 
sure (i.e. > 1) at any point along the field line. From 
the pressure, we calculate the density assuming an ideal 
gas and determine the morphology of the optically thin 
X-ray emission by integrating along lines of sight through 
the corona. We show results for two stellar inclinations, 
60° (appropriate for AB Dor) and the value of 90° deter- 
mined by [Marino et al. (2003)| for VXR45. The rotational 
modulation is defined as the fractional change in emission 
measure over one rotation. 



3. Results and conclusions 

We show in Fig. |2] images of the X-ray emission for two 
sample rotation periods, 0.514 and 0.17 days. As the scale 
height A has increased with rotation rate from 0.8i?* to 
1.3i?*, more field lines have been forced open by the pres- 
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sure of the coronal gas and there is a clear reduction in 
the volume of the corona that is bright. To quantify this, 
we define a filling factor / based on the emission-measure- 
weighted density = J rv^dV/ J n^dV such that 



JnldV 



(3) 



If the rotation period decreases from 0.514 to 0.17 days, 
/ falls from 4.7 x 10"^ to 1.8 x lO'^ while the heights of 
the brightest regions fall from about A/2 to A/10. 

As shown in Figs.|31and^there is a corresponding drop 
in the magnitude of the emission measure but a rise in its 
rotational modulation. The variation of the emission mea- 
sure for different scalings of the magnetic field strength 
and the temperature with rotation rate has already been 



explored in Jardine & Unruh (1999) The overall level of 
emission depends on where ric oz p/T oc B"^ /T, but the 
rotational modulation depends on the pattern of bright 
and dark regions. Even if the field strength and tempera- 
ture are independent of rotation rate (as in Figs. 13 and 0)) 
the rise in the pressure scale height with rotation rate will 
force closed loops to become open when 2fip > and so 
convert bright regions to dark regions. Changing the way 
that the field strength scales with rotation rate alters the 
overall level of emission, but has little effect on the rota- 
tional modulation since we have already set p cx B^ at the 
coronal base and so the gas and magnetic pressures rise to- 
gether. Lowering the temperature but keeping it and the 
field strength independent of rotation rate simply shifts 
the curves in Figs. 01 and 01 to higher values. 

The rotational modulation is less for the lower inclina- 
tion case since the brightest regions are at high latitude 
where they remain in view as the star rotates. In each 
hemisphere the brightest region is a high-latitude section 
of the large magnetic arcade that runs north-south over 
the rotation pole. Consequently, with an inclination of 90° 
these two bright regions give two maxima in the variation 
of the emission measure with rotational phase (unless they 
have been positioned at the same longitude when the com- 
posite surface map is made in which case only one peak is 
seen). For lower inclinations, the bright region in the lower 
hemisphere is hidden, and only a single peak remains. 

For the case of a 90° inclination and rotation period 
of 0.223 days appropriate for VXR45, the rotational mod- 
ulation is 26% at lO'^K and 36% at lO^K, both close to 
the observed value of w 30%. The similarity between the 



two temperatures (also noted by Marino et al. (2003) I is 
because the same structures are emitting at both temper- 
atures. Since the hotter gas has a greater scale height, 
however, it suffers more from coronal stripping and so the 
emission measure at lO'^K is about 3% of its value at lO^K. 

This model is of course not intended to reproduce ex- 
actly the observed results for VXR45. Its aim is simply to 
demonstrate the principle that rapid rotation alone can 
account for the high rotational modulation and reduced 
emission measure observed for a very rapidly-rotating star. 
At high rotation rates the increased plasma pressure in the 



summits of the tallest magnetic loops forces these loops to 
be open, with the result that X-ray bright loops become 
X-ray dark open field regions. Thus, at the highest rota- 
tion rates, much of the coronal volume is dark in X-rays 
leading to an increased rotational modulation. 
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